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ABSTRACT 

Context. We present a spectral analysis of the binary G 224-58 AB that consists of the coolest M extreme subd¬ 
warf (esdM5.5) and a brighter primary (esdK5). This binary may serve as a benchmark for metallicity measurement 
calibrations and as a test-bed for atmospheric and evolutionary models for esdM objects. 

Aims. We perform the analysis of optical and infrared spectra of both components to determine their parameters. 
Methods. We determine abundances primarily using high resolution optical spectra of the primary. Other parameters 
were determined from the fits of synthetic spectra computed with these abundances to the observed spectra from 0.4 
to 2.5 microns for both components. 

Results. We determine T e g =4625 ± 100 K, log g = 4.5 ± 0.5 for the A component and T e s = 3200 ± 100 K, log g = 5.0 ± 
0.5, for the B component. We obtained abundances of [Mg/H] = —1.51±0.08, [Ca/H] = —1.39±0.03, [Ti/H] = —1.37±0.03 
for alpha group elements and [CrH] = —1.88±0.07, [Mn/H]=-1.96±0.06, [Fe/H] = —1.92±0.02, [Ni/H] = —1.81±0.05 and 
[Ba/H] = —1.87±0.11 for iron group elements from fits to the spectral lines observed in the optical and infrared spectral 
regions of the primary. We find consistent abundances with fits to the secondary albeit at lower signal-to-noise. 
Conclusions. Abundances of different elements in G 224-58 A and G 224-58 B atmospheres cannot be described by one 
metallicity parameter. The offset of ~ 0.4 dex between the abundances derived from alpha element and iron group 
elements corresponds with our expectation for metal-deficient stars. We thus clarify that some indices used to date to 
measure metallicities for establishing esdM stars based on CaH, MgH and TiO band system strength ratios in the optical 
and H 2 O in the infrared relate to abundances of alpha-element group rather than to iron peak elements. For metal 
deficient M dwarfs with [Fe/H] < -1.0, this provides a ready explanation for apparently inconsistent "metallicities" 
derived using different methods. 

Key words, stars: abundances - stars: atmospheres - stars: individual (G 224-58 A ) - srars: Population II - stars: late 
type 


1. Introduction 


Very low-mass stars or M dwarfs are the most numerous and 
longest-lived stars in our Milky Way. They are targeted for 
the study of the Galactic structure and searches of habitable 
exoplanets by a large number of projects. 

The precise metallicity meas urement of M dwa r fs has 
become a popular field (e.g., iRoias-Avala et all 120101 : 

f Onehag et al.l 2012 : Neves et~aLll2013l; 
However, the data quality of obser- 
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vations and models make it useful to analyse M dwarfs in 
binary systems with F-, G- or K-type dwarfs in order to en¬ 
sure robust calibration of metallicity features in their spec¬ 
tra. These M dwarfs not only share the same age and dis¬ 
tance with their FGK companions. They also share chem¬ 
ical abundance which have been well calibrated for FGK 
dwarfs. Such studies are limited to the M dwarfs in the 
Galactic disk where there exist enough wide FGK+M dwarf 
binary systems for calibration. 


Send offprint requests to: Ya. V. Pavlenko 


Sub classi fications for low metallicity M dwarfs or 
M subdwarfs dGizisI 119971 : iLepine et al.l I2007T) are widely 
in use. Precise metallicity measu rements of M subd warfs 
have been recent l y con ducted by IWoolf et al.1 ((2009) and 
iRaipurohit et all (l2014fl . They both used high resolution 
spectra of subdwarfs to measure metallicities and calibrate 
a relation between [Fe/H] and molecular band strength in¬ 
dices from low-resolution spectra. This calibration makes 
possible to estimate the metallicity of a large sample of 
subdwarfs without the observational cost (some times im¬ 
possible) of high resolution spectroscopy of th ese faint ob¬ 
j ects, to enable meaningful population analysis. IWoolf et al.l 
(j2009f l measured [Fe/H] from Fe I lines using NextGen 
models to fit the spectra of 12 low-metallicity main se¬ 
quence M stars. Adding these objects to previous c ompi- 
lations (IWoolf fe Wallersteinl 120061 : IWoolf et al.ll2009D they 
established the relationship between the CTiO/CaH param¬ 
eter and metallicity for the objects with T P f f o f 3500 -4000 
K, and metallicities over ^-2. IRaipurohit et al.l (12014D used 
BT-Settl models to fit the spectra of three late-K subdwarfs 
and 18 M subdwarfs (covering the M0-M9.5 range), reach- 
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ing the lowest value of [Fe/H]=-1.7. These works mainly 
use the spectra of single M subdwarfs, without the support 
of the measuremen ts of any binary FGK comp anion. It is 
worth noting, that IWoolf fc Wallerstein (l2006f) studied M 
dwarfs with a G or warm K companion and M dwarfs in 
the Hyades cluste r as part of their sample. Furthermore, 
iBean et~akl ( 2006 1 reported tests of metallicity measure¬ 
ments using M dwarfs in binaries with warmer companions. 

G 224-58 AB is a esdK5+esdM5.5 wide binary system 
with 93 a rcsec separati on identified by motion, colour and 
spectra in IZhang et al.l (|2013T) . G 224-58 B is the coolest M 
extreme subdwarf with a bright primary. In this paper, we 
study G 224-58 A through the analysis of high resolution 
observed and synthetic spectra. This provides by associa¬ 
tion a more precise metallicity of G 224-58 B a very cool 
extreme subdwarf that will serve as benchmark for cali¬ 
bration of metallicity measurement and as a test-bed for 
atmospheric and evolutionary models of such metal-poor 
low-mass objects. In Section [5] we present the observational 
data used in our analysis. Procedure and results of our anal¬ 
ysis are given in Sections [3] and [U respectively. Our main 
results are discussed in section [5] 

2. Observations 

2.1. FIES optical spectrum of G 224-58 A 

High resolution spectrum of G 224-58 A were taken on 20 13 
June 19 with the the high-resolution (FIES; iTelting et alJ 
120141) 1 on the Nordic Optical Telescope. The CCD13, 
EEV42-40 detector (2kx2k) was used in medium resolu¬ 
tion mode (R=46000) covering the 3630 - 7260 A interval 
in 80 orders. Signal-to-noise is ~35 at Ha. In general, this 
value is a bit low for classical abundance analysis using the 
curve-of growth procedure. In our case only well defined 
and fitted parts of spectral lines are used. Another problem 
is the determination of the continuum level. 

Data were reduced using the standard reduction proce¬ 
dures in IRAFfl We used the echelle package for bias 
subtraction, flat-field division, scatter correction, extrac¬ 
tion of the spectra, and wavelength calibration with Th-Ar 
lamps. Since FIES have a dedicated data reduction pipeline, 
FIESTooB we compared our reductions with those of the 
pipeline and found them to be equivalent. 

2.2. SDSS optical spectrum of G 224-58 B 

As in IZhang et alJ (1201311 . we used the SDSS 
J151650.33+605305.4 (G 224-58 B) spectrum from 
the Science Archive Server (SASB that contains the 
Tenth SDSS Data Release (DR10) of Sloan Digital Sky 
Survey (SDSS; lYork et all 1 2000 1. The data is provided 
as final spectrum product from the pipeline. The SDSS 
spectroscopic pipelines, extract one dimensional spectra 
from the raw exposures produced by the spectrographs, 
calibrate them in wavelength and flux, combine the 
red and blue halves of the spectra, measure features in 

1 IRAF is distributed by the National Optical Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under contract with the National Science 
Foundation. 

2 http://www.not.iac.es/instruments/fies/fiestool/FIEStool.html 

3 http://drl0.sdss3.org/ 


these spectra, measure redshifts from these features, and 
classify the objects as galaxies, stars, or quasars (see 
http://www.sdss3.org/drl0/spectro/pipeline.php for more 
details). 

2.3. LIRIS near infrared spectra of G 224-58 AB 

Medium resolution near infrared (NIR) spectra of G 224-58 
AB were obtained with the Long-slit I ntermediate Resol u- 
tion Infrared Spectrograph (LIRIS [Manchado et alJ 11998T) 
at the William Herschel Telescope on 2015 May 3 and 4. 
We used hr_j, hr_h and hr_k grisms for the observation, 
which provide a resolving power of 2500. The data were re¬ 
duced with the IRAF LIRISDR package which is supported 
by Jose Antonio Acosta Pulido at the IACB The total in¬ 
tegration times used on G 224-58 A were 600s divided into 
ten 60s exposures using the technique of ABBA exposures 
to remove sky background in J, H and K bands. Three 
of the K band spectra of G 224-58 A were taken on the 
gap between detectors and rejected for final combination. 
Thus the total integration time for final K band spectrum 
of G 224-58 A is 420s. The total integration times used on 
G 224-58 B were 18 x 120s = 3600s in J and H bands and 
20 x 120s = 4000s in K band. 


3. G 224-58 A spectrum analysis 

3.1. The procedure of the G 224-58 A abundance analysis. 

To determine the physical parame ters of G 224-5 8 A w e 
followed the procedure described in lPavlenko et all (|2012f ). 
We determined abundances in the G 224-58 A atmosphere 
iteratively. Namely, after the determination of abundances, 
the model atmosphere was recomputed to account for the 
changes to these abundances. The ID model atmosp heres 
were computed by SAM12 program ( Pavlenkol 12003 1. Ab¬ 
sorption of atomic and molecular lines was accounted for 
by opacity sampling. Spectroscopic data for absorption 
lines of atoms and molecules were taken fr om the VALD 2 
dKupka et al.l Il999h and Kurucz database (iKurucdl19931) . 
respectively. The shape of each atomic line is determined 
using a Voigt function and all damping constants are taken 
from line databases, or computed using Unsold’s approach 
(Unsol d 1954). For th e complete description of routines (see 
iPavlenko et al.ll2012h . 

Sy nthetic spectra are calculated with the WITA6 pro¬ 
gram (|Pavlenkolll997 ), using the same approximations and 
opacities as SAM12. To compute the synthetic spectra we 
used the same line lists as for the model atmosphere compu¬ 
tation. A Wavelength step of AA = 0.025 A is employed in 
the synthetic spectra computations. The synthetic spectra 
that are computed across the selected spectral region are 
convolved with profiles that match the instrumental broad¬ 
ening and that take into account r otation al broadening fol¬ 
lowing the procedure described bv iGravl ( 1976il . 

We used the Sun as a template star in this study ap¬ 
plying our technique to determine solar abundances in or¬ 
der to give reliability to our abundances for G 224-58 A . 
We used an observed solar spect rum with resolution R > 
100,000 from lKurucz et all (Il984i) . The abundances are de¬ 
termined for a list of pre-selected lines which are fitted well 

4 http://www.iac.es/gaIeria/jap/lirisdr/LIRIS_DATA_RE- 
DUCTION.html 
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Fig. 1 . Top: fits of the observed in spectrum G 224-58 A Mg 
I line triplet and MgH lines Labels ’atoms’ and ’atoms+MgH’ 
mark synthetic spectra computed with account of atomic lines 
and atomic lines + MgH, respectively. Middle: dependence of 
Fk on Vrn for Fe lines, model atmosphere 4625/4.5. Bottom: fits 
to Na I resonance doublet. 


in the solar spectrum and for the Sun we adopt T 
g/[Fe/H] =5777/4.44/0.00 (see iPavlenko et al.ll20ll V The 
abundances we obtained are in good agreeme n t wit h 
known classical resu l ts, e. g. lAnders fe Grevessel (11989H ; 
iGurtovenko fe Kostikl (|l989f) within ± 0.1 dex. Further de¬ 
tails of the procedure are given in lKushniruk et all (l2014lf . 


3.2. Modelling G 224-58 A optical spectrum 

From the literature, we expect that the effective temper¬ 
ature of G 224-58 A is in the 4500 - 4750 K range and 
that the atmosphere is metal deficient. We carried out our 
abundance determination of the star applying the proce¬ 
dure described in the section EEJ for model atmospheres of 
T e ff = 4500 and 4750 K, log g = 4.5 and 5.0. These compu¬ 
tations provide rather similar results, i.e [Fe/H] ~ -2.0. We 
probed these solutions with analysis of the spectral region 
where strong lines of Mg I and the A 2 II — X 2 E + band 
system of MgH molecule are located, see the upper panel of 
Fig. [TJ These features in the synthetic sp ectra show a no¬ 
table dependence on the gravity chosen (|Kushniruk et al.1 
120131) . Additionally we determined the titanium abundance 
log N(Ti) using the lines of different titanium ions, i.e. Ti I 
and Ti II. We found formally similar results of abundance 
determinations for the cases of 4500/4.5 and 4750/5.0. How¬ 
ever, MgH lines were too strong in the case of 4500/4.5 and 
too weak in the case of 4750/5.0, and abundances obtained 
from the fits to Ti I and Ti II lines disagree by more than 
0.1 dex with each other, see Section Hi. 2. 31 

We checked T e g = 4625 K, and log g = 4.5-5.0 mod¬ 
els. Since 4625/4.5 model matches better both the MgH 
lines and the minimal disagreement between Ti I and Ti II 
abundances (see section l3.2.3p . we choose it for our further 
analysis. 

3.2.1. Microturbulent velocity V m 

The procedure of microturbu lent velocity determ ination 
differs from that described in IPavlenko et ali ( 2012 1 as fol¬ 
lows. 

— We determine a small grid of N v microturbulent veloc¬ 
ities Vm = Vo + i * 0.5, i changes from 1 to 7. We adopt 
Vo = 0 km/s and so cover a grid of plausible values for 
V m in the atmospheres. 

— For every microturbulent velocity from our grid we ob¬ 
tained the iron abundances logNk{Fe), k = 1,... L, here 
L is the number of lines in our pre-selected list, and de¬ 
termined the central residual fluxes of fitted lines r^k. 

A o 

— Using x 2 fitting, we approximate the dependence of log 
N(Fe) vs. T\ 0 by the linear formula logN(Fe) = A + 
F a * (r\\ 0 ), here A , F a are constants which depend on 
Vm . We are interested in the case F a = 0, when all lines 
provide the same abundances. 

A similar approach for V m determination from the 
equivalent width vs. abundance analysis was used by R. Ku- 
rucz in his WIDTH9 code. However, our approach is more 
flexible in the sense that we can use V m determination even 
for blended lines. 

In the middle panel of Fig. |T| we see F a = 0 at least 
at V m =0.5 and 2.5 km/s. Corresponding abundances of 
iron obtained from the fits to profiles of 44 lines of Fe I in 
the spectrum of the primary are -6.290 ± 0.021 and -6.478 
± 0.026, respectively. We choose the first value due to the 
lower dispersion of results. It is worth noting, that V m = 
0.5 km/s seems to be reasonable from the point of view 
that the lower opacity photospheres of metal poor stars al¬ 
low spectra to probe deeper, i.e. into higher pressure layers 
where we cannot expect high velocity motions/large V m . 
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Table 1. Fitted Ti II lines in spectra of G 224-58 A and the 
Sun, spectroscopic data were taken from VALD3. 


Wavelength, A 

gf 

E"(eU) 

5129.16 

4.571E-02 

1.892 

5226.54 

5.495E-02 

1.566 

5336.79 

2.512E-02 

1.582 

5381.02 

1.072E-02 

1.566 

5418.77 

7.413E-03 

1.582 


3.2.2. Fit to sodium resonance doublet 

Sodium lines are of special interest as they form strong 
features in the spectra of both G 224-58 A and G 224- 
58 B . In the primary spectrum we see the strong lines 
of the resonance doublet of Na I at 5889.95 and 5895.92 A. 
Their intensities depend on T e g , log g , V m and the sodium 
abundance. In the bottom panel of Fig |T| we show the fit 
of our synthetic spectra to the Na I doublet line profiles 
observed in G 224-58 A spectrum. Both lines are better 
fitted with the model atmosphere T e g = 4625, log g =4.5 
with N(Na) = -7.34, i.e, we obtain [Na/Fe] = -0.2 for the 
primary. 


3.2.3. Fit to Ti II lines in the observed spectrum of 
G 224-58 A 

The abundance ratio of ionized species of the same element, 
usually Fe I/Fe II, is known to change with log g for a fixed 
T e g and therefore it is often used to verify the selection of 
log g , serving to constrain the global solution. However, 
our primary G 224-58 A is metal deficient and a relatively 
cool star. In comparison with the Sun, it is difficult to locate 
appropriate lines of Fe II. We found a few features consist¬ 
ing of weak and blended Fe II lines, but results of their fit 
was not confident due to the low quality of the observed 
spectrum. 

On the other hand, due to the lower ionization poten¬ 
tial of the neutral titanium (E = 6.83eU) in comparison 
to Fe I (E = 7.90eU), the lines of Ti II should be stronger 
in spectra of comparatively cooler stars, i.e., more useful 
for abundance analysis. Therefore we used lines of Ti I and 
II to make the analysis instead. Indeed, we found lines of 
Ti II which are fit well in our synthetic spectra and list in 
Table [1] To be more confident, we also looked at these lines 
in spectrum of the Sun finding a good agreemen t between 
the lin es in model atmosphere 57 77/4.44/0 fromlPavlenkol 
(|2003lf and observations from the iKurucz et al.1 (11984f) at- 
las, see Fig. [2j The good agreement between abundances 
found using Ti I vs Ti II gives additional independent sup¬ 
port in our choice for the G 224-58 A model atmosphere 
parameters. 


3.2.4. Other elements 

Abundances determined from selected line lists are shown 
in the Tabled In the 6-th column of the table we show the 
number of preselected line lists and number of accounted 
lines. In particular we note that that light a—elements Mg, 
Ca, Ti are overabundant in comparison with elements of 
iron peak Cr, Mn, Fe, Ni, Ba. For this analysis we use the 


weaker lines, i.e not ones that we find to be saturated in the 
observed spectrum. Table [5] indicates that the uncertainties 
of the abundances do not exceed 0.17 dex for the adopted 
range of parameter variation. 

The found overabundance of Mg, Ca, Ti with respect to 
the Fe shows rather weak dependence on the choice of T e g , 
log g or V m , see Table [5] In general, the overabundance 
of alpha-elements is a well kn o wn p he nomenon fo r meta l- 
deficien t star s, e.g., iMagainl (I1987H: iMcWilliaml Il997fl : 
ISnedenl (I2004T) . or more recentlv lHansen et al.l (I2015T) . 

To carry out the process of abundance analysis we fitted 
the observed profiles of the spectral lines. This procedure 
allows us to determine the rotational velocity v sin i . It is 
worth noting that the values of v sin i in Table 2 should 
be considered rather as upper limits due to the natural re¬ 
strictions provided by the limited quality of our observed 
spectrum. In some sense v sin i is used here as the adjust¬ 
ing parameter of our fitting procedure. For more a accurate 
determination of v sin i we should adopt/develop more so¬ 
phisticated models of instrumental broadening and macro¬ 
turbulence and use them for fits to the observed spectra of 
better quality. For now we can claim that G 224-58 A is 
a slowly rotating (v sin i < 3 km/s) star, this agrees well 
with its status of old halo star. 


3.3. LIRIS infrared spectrum of G 224-58 A 

Fits to LIRIS low resolution (R=2500) spectra observed 
in H , J, K spectral ranges are shown in Fig. [3] The Fig¬ 
ure shows the observed spectra, where we overplotted the 
theoretical spectra (blue and green lines), for two different 
Silicon abundances. The theoretical spectra were computed 
by using the optical G 224-58 A spectrum (see section 15721) . 

Despite the comparatively low resolution, we are able to 
identify a few Si I lines in the observed spectra. In the Fig. [3] 
we show two cases, with [Si] =-1.0 and -2.0. It is apparent 
that Si=—1.0 is a better fit and thus Si I is overabundant. Si 
(Z=14) belongs to the alpha-element group and as discussed 
earlier its overabundance is expected. Here we can claim 
only an estimation of Si overabundance, i.e. [Si/H] ~ -1.1 
±0.3 dex based on an average of the line fits. Here we note 
that changes of Si abundance by factor 10 affects the general 
shape of the observed SEDs (Fig [3]). Si I is an important 
donor of free electrons in cool atmospheres, therefore both 
model atmospheres and computed spectra show dependence 
on its abundance. Unfortunately, the absorption lines of Si 
I are too weak to discern in our optical spectrum of G 224- 
58 A lower limit. 

In the H spectral region we see a strong Mg I line at 
15754 A. We obtained a fit to this line with the Mg abun¬ 
dance obtained from the analysis of optical spectrum. In 
fact, this provides independent confirmation of our results 
in Table 2. 

In the K spectral region we can discern the presence 
of the CO overtone bands at the expected level but not 
sufficiently well to determine the carbon abundance. 


4. G 224-58 B spectrum analysis 

4.1. Procedure 

G 224-58 B is a late spectral type dwarf. In general we will 
assume that, as a binary the B component has the same 
properties of age and composition that we have found for 
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Fig. 2. Fit of our theoretical Ti II lines to the observed spectra of the G 224-58 A and the Sun. In the theoretical spectra only 
absorption by Ti II is included. Spectra of G 224-58 A are shifted vertically to simplify the plots. 


Table 2. Abundances of elements determined from fits to G 224-58 A hires spectrum. I and L to tai are a number of accounted 
lines and the total number of lines in pre-selected line list L, respectively. Sensitivity of abundances on the changes of parameters 
of model atmosphere with respect to the model 4625/4.5/-1.5*, V m = 0.5 km/s model is shown in the latter columns. 


Elem 


log N(x) 

|X/H| 

V 

sin i 

l / Ltotai 

A V m = +0.5 km/s 

A log g = 0.5 

A T eff = -100 K 

Mg 

12 

-5.972 

-1.512 

2.33 

± 

0.42 

6/7 

0.003 

-0.007 

0.087 

Ca 

20 

-7.067 

-1.387 

2.66 

± 

0.14 

16/16 

0.038 

0.012 

0.107 

Ti 

22 

-8.422 

-1.372 

2.25 

± 

0.21 

14/16 

0.029 

-0.068 

0.114 

Cr 

24 

-8.254 

-1.884 

2.83 

± 

0.12 

9/9 

0.089 

0.033 

0.167 

Mn 

25 

-8.609 

-1.959 

2.50 

± 

0.32 

5/5 

0.040 

-0.100 

0.080 

Fe 

26 

-6.290 

-1.920 

2.50 

± 

0.10 

44/44 

0.051 

-0.094 

0.090 

Ni 

28 

-7.598 

-1.808 

1.89 

± 

0.32 

9/9 

0.089 

-0.134 

0.092 

Ba 

56 

-11.780 

-1.870 

1.88 

± 

0.38 

4/4 

0.125 

-0.125 

0.166 


A component. However, we perform a different approach to 
metallicity determination. 

First, we generated a grid of the LT E synthet i c spec - 
tra following the procedure outlined in IPavlenko (t 1997T ). 
We started with the NextGen model atmospheres of ef¬ 
fective temperatures 2800, 2900, 3000, 3100, 3200, 3300, 
3400 K, gravities of log g = 4.0, 4.5, 5.0, 5.5 and metal- 
licities [Fe/H] = -1.5 and 2.0. We replaced some abun¬ 
dances by the values obtained for the G 224-58 A with 
parameters T e g =4625 K, log g =4.5. In our computa¬ 
tions the line lists of VO and CaH were taken from Ku- 
rucz’s websit(0, for more details see iPavlenkol (120 14 ). The 
CrH a nd FeH line lists were c omputed by iBurrows et al.l 
(|2002fl and iDulick et all (l2003lf. re spectively. We upgraded 
the TiO line lists of Pled ( 19981) with the new version 
available on his websitcjfl- Infrared spectra were computed 
with account of water v apour absorp tion lines provided 
by the EXOMOL group iBarber et all (120061) . The spec¬ 
troscopic data for atomic absorptio n VALD3 c o me fro m 
the Vienna Atomic Line Database iKuoka et al.l (Il999lf1 . 
The profiles of the Nal and Kl resonance doublets were 
computed here in t he framework of a q uasi-static ap¬ 
proach described in [Pavlenko^t_alJ_(200^) with an up¬ 
graded approach from IBurrows fe Volobuvevl (|2003l) . More 


5 kurucz.harvard.edu 

6 http://www.pages-perso-bertrand-plez.univ-montp2.fr/ 

7 http://vald.astro.univie.ac.at/ vald/php/vald.php 


details on the techni que and procedure are presented in 
IPavlenko et al.l (120061) . 

4.2. Results of the G 224-58 B optical spectrum analysis 

In Fig. 2] we show a compilation of the main absorption fea¬ 
tures that can be observed in the spectrum of G 224-58 B . 
The main feature is formed by CaH A 2 H r — V 2 E + band 
system. Strong MgH, TiO band systems and K I and Na 
I resonance lines also form noticeable features in the sec¬ 
ondary spectrum. 

The choice of best fit between the observed and com¬ 
puted spectra was achieved by minimizing the parameter 

s{h) = F ?) 2 ’ 

V 

where F v and F* are the fluxes in the observed and com¬ 
puted spectra, respectively, and // i s the normalizat i on fac - 
tor. A similar procedure was used bv IPavlenko et ab (120061) . 
Fits were made for all synthetic spectra in our grid. We use 
the instrumental resolution R = 1500 for our calculations. 

The best fit of our synthetic spectrum to the G 224- 
58 B observed SED is obtained for the 3200/5.0/-1.5* model 
atmosphere (see Fig. [5]) - Abundance scaling parameter -1.5* 
means here that a reducing factor -1.5 was used for all abun¬ 
dances in the NextGen model atmosphere except for those 
obtained from the analysis of the primary spectrum. Then, 
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Wavelength (A) 


Fig. 3. Fit to G 224-58 A LIRIS infrared spectrum (top left). Three others plots show fits to the observed H (top right), J (bottom 
left), K (bottom right) spectral regions in larger scale 



Fig. 4. Left: identification of the main molec ular features formed th e observed optical spectrum of the B component, model 
atmosphere 3200/5.0/-1.5* from NextGen grid (lHauschildt et ah (1999)). Right: fit of our synthetic spectrum computed with the 
same model atmosphere to the observed SED of G 224-58 B . 


we used ’’the best fit” (as defined previously) to describe the 
best agreement between the computed and observed SED 
and the profiles of strong lines. Fig. [5] shows the fit to the 
observed SEDs and intensities of the strong sodium and 
potassium lines, with abundances obtained from analysis 
of G 224-58 A spectrum. 

4.2.1. Strong features in the G 224-58 B spectrum 

The notable features in the observed spectrum are formed 
by MgH A 2 II — X 2 £ + band system at 5200 A, CaH 


A 2 n r — Af 2 E + and B/B'2Y, + — X 2 E + band systems in the 
wide spectral range, as well as TiO 7 band system at 7200 
A. We see that these features in the observed spectrum are 
well reproduced by our theoretical computations with the 
abundance obtained from G 224-58 A analysis. 

Furthermore, in the observed G 224-58 B spectrum we 
see sodium and potassium resonance doublet lines at 5890 
and 7680 A, as well as sodium triplet lines at 8190 A known 
as a good gravity discriminator. In Fig[5]we show our fits to 
these observed K I and Na I lines. To reproduce Na I lines 
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Wavelength (A) Wavelength (A) 




Wavelength (A) Wavelength (A) 

Fig. 5. Fits of some spectral details in our theoretical spectra shown in the right panel of the Fig 0] to the observed features in 
G 224-58 B spectrum: MgH A 2 II — X 2 S + band system at 5200 A(top left), lines of Na I resonance doublet (top right), lines of 
resonance doublet of K I (bottom left), subordinate lines of Na I (bottom right). Model atmosphere and normalization point are 
the same as for Fig [4] 


in G 224-58 B spectrum we used the sodium abundances 
determined from analysis of G 224-58 A spectrum. 

The potassium resonance lines are beyond the spectral 
range observed for the G 224-58 A , and other potassium 
lines are very weak. Still we see strong lines of potassium 
doublet in G 224-58 B observed spectrum. To fit these lines 
we used the same abundance scaling factor, as was found 
for sodium, i.e. [K/Fe] = -0.2 and obtained good agreement 
between observed and computed profiles. It would be inter¬ 
esting to verify this result by the comparison of computed 
and observed profiles of potassium resonance doublet in the 
G 224-58 A spectrum. 


4.2.2. Fe lines 


Since iron lines in G 224-58 A spectrum are numerous, the 
determination of Fe abundance was done with high accu¬ 
racy. Ideally a similar analysis would be performed for the 
G 224-58 B , using the same Fel line list. However, practi¬ 
cally its lower effective temperature, where molecular bands 
dominate the spectrum combined with the lower resolution 
and signal-to-noise make this difficult. The lines are at the 
level of the noise. We could use only the Fe I line at 5325.19 
A. We only can claim that weak Fe I lines are present in 
the observed spectrum and based on Fig 0 find [Fe/H] = 
-1.7. 


4.3. Fit to LIRIS spectra of G 224-58 B 

The LIRIS infrared spectrum of G 224-58 B has compara¬ 
tively low signal-to-noise. Therefore, to reduce the level of 
noise we binned the observed spectrum by a factor of 7. In 
the following we show the binned spectra to simplify. 

Synthetic spectra were computed for the model atmo¬ 
sphere 3200/5.0/-1.5* with a 0.05 A step taking into ac¬ 
count all known molecular opacity sources. First of all we 
computed the contribution of the molecular bands to the 
opacity in our spectral ranges. Results are shown in the left 
panel of Fig [5] Water vapour absorption dominates across 
all modelled spectral ranges. In the most general case ab¬ 
sorption by H 2 O depends on effective temperature, oxygen 
abundance and other input parameters. In our work we 
fixed almost all of these parameters, so we have a chance to 
estimate the oxygen abundance. Indeed, the shape of the 
computed SEDs across the H band depends on the water 
absorption, with better agreement when we adopt some un¬ 
derabundance in the atmosphere of the B component, see 
the right panel of Fig. [HI 

4.4. Fit of BT-Settl spectra to the optical and infrared SEDs 
of G 224-58 B 

Here we computed synthe tic spectra for model at mospheres 
from the NextGen grid (lHauschildt et all I1999D as newer 
BT-Settl model atmosphere are not in the public access. 
Nonetheless, fluxes for a fixed grid of abundances are avail- 
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Fig. 6. Left: identification of the main molecular features formed the observed by LIRIS infrared spectrum of the B component, 
model atmosphere 3200/5.0/-1.5*. ”Vald” marks the synthetic spectrum computed for atomic line list. ”H”, ”J”, ”K” mark three 
spectral ranges observed by LIRIS. Right: fit of our synthetic spectra computed with the same model atmosphere, but different 
oxygen abundances, to the observed SED of G 224-58 B . 



Wavelength (A) 


Fig. 7. Fit of our computed 3200/5.0/-1.5* model atmosphere 
theoretical spectrum to the observed G 224-58 B spectrum. NIR 
spectra were smoothed here by 3 pixels for clarify. Vertical ar¬ 
rows mark positions of Fe I lines in the computed spectrum. 


able and we used them for the comparison with our ob¬ 
served spectra. 

In Fig. [ 8 ] we provide fits of the BT-Settl spectra com¬ 
puted with 3150/5.0/-1.5 model atmosphere to the observed 
fluxes of G 224-58 B . Here the metallicity parameter [Fe/H] 
= -1.5 was used to adjust all abundances. As can be seen, 
BT-Settl fluxes reproduced well the observed SEDs. Also 
the Ca, Mg and Ti abundances adopted in BT-Settl model 
agree well with our values, obtained from the fits to the A 
component. Therefore CaH, MgH and TiO bands are fit¬ 
ted well enough with both our model with modified abun¬ 
dances and an original BT-Settl model. We should remark 
that the NextGen and BT-Settl model atmospheres should 
be similar. Indeed, dusty effects at these relatively warm 
temperatures from the BT-Settl model atmosphere will be 
very minor, given the case of similarity of the main opacity 
sources we should get similar results. 


Table 3. Abundances determined for the atmospheres of G 224- 
58 A and G 224-58 B . The errors given for G 224-58 A are de¬ 
termined by our minimisation procedure and for G 224-58 B are 
made “by-eye”. 


Element 

Z 

log N 

[x/h\ a 

[X/H\ b 

O 

8 

-5.3 




- 2.1 

± 

0.2 

Na 

11 

-7.87 

- 2.1 

± 

0.2 

- 2.1 

± 

0.2 

Mg 

12 

-5.975 

-1.512 

± 

0.079 

-1.5 

± 

0.1 

Si 

14 

-5.59 




- 1.1 

± 

0.3 

K 

19 

-9.0 




- 2 . 1 ± 

0.2 

Ca 

20 

-7.105 

-1.387 

±0.025 

-1.4 

± 

0.2 

Ti 

22 

-8.451 

-1.372 

± 

0.028 




Cr 

24 

-8.343 

-1.884 

± 

0.065 




Mn 

25 

-8.649 

-1.959 

± 

0.055 




Fe 

26 

-6.341 

-1.920 

± 

0.020 

-1.9 

± 

0.5 

Ni 

28 

-7.687 

-1.808 

± 

0.050 




Ba 

56 

-11.905 

-1.870 

± 

0.111 





5. Discussions and conclusions 

Our final results of abundance determination in atmo¬ 
spheres of G 224-58 A and G 224-58 B are shown in Table 
[3] Here we compiled abundances obtained by fits to our 
spectra across optical and infrared spectral ranges, using 
different approaches to get the best fits to observed absorp¬ 
tion profiles and spectral features. In that way we obtained 
the fits to G 224-58 A and G 224-58 B spectra using the 
same abundances. 

IZhang et al.l (|2013l ) provided very strong arguments 
about the binarity of G 224-58 AB, based on the common 
high proper motion (PM) and radial velocity (RV). We used 
their results as a starting point for our investigation. Our 
study provides confirmation about the binary nature of the 
objects through the consistency of the abundance measure¬ 
ments of the A and B components. 

From our analysis of the primary spectrum (see Ta¬ 
ble 0 we determined [Mg/Fe] = +0.41±0.05, [Ca/Fe] = 
-f0.53±0.05, [Ti/Fe] = +0.55±0.05, i.e. these elements are 
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Wavelength (//m) 


Fig. 8. Spectra of G 224-58 AB compared to BT-Settl model spectra. Spectra are normalised at 0.82 pm. 


overabundant relative to iron. Our finding is confirmed 
by direct modeling of the cool secondary dwarf spectra. 
Generally speaking, the overabundance of Ca in atmo¬ 
spheres of metal poor stars is w ell known (see iMagainll 1 98?t 
lMcWilliamlll997t ISnedenll2004 and references therein). In 
the case of binary systems, or exoplanetary systems the sit¬ 
uation looks even more intriguing, because here we have the 
impact of a few very different processes: early epoch nucle¬ 
osynthesis, separation of heavy elements in the binary sys¬ 
tem process formation, etc. An interesting aspect of study¬ 
ing halo systems including a cool component like G 224- 
58 AB is that it provides an opportunity for the abundance 
determination of some elements, e.g. K and Na. 

Generally speaking, for [Fe/H] > -1.0 there is not much 
variation in abundance ratios for different elements at any 
given [Fe/H] for single (non-binary) main sequence stars, 
so that knowing [Fe/H] does give us a good idea what its 
ab undances are f or oth er elements, e.g., figures 9, 10, and 11 
in iReddv et al.1 (l2003i ). Besides, for lower metallicity stars 
there is more scatter in the [X/Fe] ratio for some elements 
for a given [Fe/H], which may affect the stellar atmosphere, 
e.g. two stars with [Fe/H] = -1.5 but with very different 
alpha element abundances will have different line opacity. 

The essential result of our G 224-58 A spectrum analysis 
is the conclusion that at least for the case of metal deficient 


stars with [Fe/H] < -1.0 we cannot use one parameter, the 
metallicity, to describe the behaviour of individua l abun¬ 
dances. This conclusion agrees well with lWoolf et al.l (I2009T) 
results and demonstrates the importance and necessity of 
abundance analysis of binaries li ke G 224-58 AB. Compar¬ 
ison of our results for Ti with I Woolf et al.l (|2009f ) abun¬ 
dances for esdM shows at least a good qualitative agree¬ 
ment. For the most metal poor esdM LP 251-35 (3580/5.0/- 
1.96) they obtained [Ti/Fe] = 0.45. Interestingly, these au¬ 
thors used K I and Ca II lines to determine the parameters 
of the atmosphere of the stars. Here we have obtained an 
overabundance in calcium and under abundance of potas¬ 
sium (and sodium) in the atmosphere of G 224-58. If the 
overall abundances of the extremely metal-poor M-dwarfs 
are similar, their results may be affected by abundance un¬ 
certainties. 

iLepine et al.l ( 200 7) defined a metallicity index, 
CTiO/CaH; which is formed by a combination of the 
spectral indices of Ti05, CaH2, and CaH3. They clas¬ 
sified M subdwarfs into three metal classes: subdwarf 
(sdM; 0.5 < OriO/CaH < 0.825), extreme subdwarf 
(esdM; 0.2 < CTiO/CaH < 0-5 ) and ultra subdwarf 
(usdM; CTiO/CaH < 0.2). Metallicity measurements of 
M subdwarf s have been conducted with high resolu- 
tion spectra dWoolf fc Wallersteinl 12006 :: IWoolf et al.ll2009 : 
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iRaipurohit et all I2014I) . The relation ship between the 
CTiO/CaH and [Fe/H] has been studied (I Woolf et al.l 120091 
iMann et al.l 120131 : IRaipurohit et al.ll2014f) . We calculated 
the metallicity index CxiO/CaH = 0.22 for G 224-58 B to 
compare to previous works. 

Figure [U] shows the CxiO/CaH and [Fe/H] of G 224- 
58 B and M subdwarfs from IWoolf et al.l i 2 (1(1!) ). A blue 
dashed line in Figure [HI represents the correlation of 
Cxio/CaH an d [Fe/H] for late-type K and early-type M sub¬ 
dwarfs (3500-4000 K) ([Woolf et al.l 12003 *1. A black solid 
line repre sents our fit t ing o f G 224-58 B and M subd¬ 
warfs from lWoolf et al.l (l2009t) . which can be described with 
[Fe/ H] = 2.004 x^xio / caH — 1-894. We excluded objects 
from IRaipurohit et al. (i2014f) in our fitting because to re¬ 
duce the dispersion of results in Fig. [9] would re quire a shift 
of [Fe/H] values from iRaipurohit et al.l (|2014| ) by a factor 
of -0.3 to -0.5 downward. This could arise from the over¬ 
abundance of Ti and Ca and their inclusion in the [Fe/H] 
values they compute. 

It is worth noting that: 

— IRaipurohit et all (]2014f l measure [Fe/H] with only iron 
lines, however, for thre e M dwarfs shown in fig. 11 of 
IRaipurohit et al.l ([2014P they find CTiO/CaH > 0.825 for 
metallicities from —0.5 to —1.0 and so these values are 
significantly discrepant for the CxiO/CaH parameter. We 
note their fit to the sdM9.5 object in a small region of 
optical spectrum gave [Fe/H] = -1.1 and T e g = 3000K. 
From the fit to the 0.7-2.5 //m NIR spectrum of the same 
object with same model, we got [Fe/H] = -1.6 and T e g = 
2700K. We note that an analysis of a narrow spectral 
region can be significantly impacted by relatively dif¬ 
ferent pseudocontinuums formed by molecular bands in 
the optical and the infrared. The measurement of abun¬ 
dance and temperature from atomic line strength can 
be significantl y impa cted by such differences, e.g., (see 
iPavlenko et ail 1 19951) . 

— Woolf’s results of CTiO/CaH determination are affected 
by uncertainties in the adopted K, Na, Ca and Mg abun¬ 
dances. 

— Our value of iron abundance obtained for the G 224- 
58 A in the analysis explained in previou s sections is 
[Fe/H ] = -1.920; only one M subdwarf from lWoolf et alJ 
(l2009h has lower [Fe/H] than G 224-58 B but was ex¬ 
cluded in their fitting. 

For the case of halo dwarfs of later spectral classes, i.e. 
esdM dwarfs, the use of one parameter of metallicity may 
seriously impact results. And, this is why we need binaries 
to calibrate spectral indices vs. abundance dependence for 
M subdwarfs of different populations. 

In general, our fits to the infrared spectra of the G 224- 
58 A and G 224-58 B observed with LIRIS provide the in¬ 
dependent confirmation of the correctness of our choice of 
parameters of atmospheres of both components. With the 
spectrum of G 224-58 B , we find evidence of overabundance 
of Si I and underabundance of oxygen. It is worth nothing 
these results were obtained by fitting to low resolution in¬ 
frared spectra. In spectra of stars of normal metallicity the 
results are affected by the presence in the spectral analysis 
of complex blending. However, in the case of metal deficient 
stars these effects are less pronounced. Moreover, in our case 
we have not any other chance to estimate oxygen and silicon 
abundances, except through the modelling of infrared spec¬ 
tra. To analyse weak oxygen and silicon lines in the optical 



Cl5 O/CaH 

Fig. 9. Metallicity parameter Cx io/PaH vs. |Fe/H| for G 224 - 
58 B and M subdwarfs from IWoolf fc Wallersteinl (|2006l ). 

I Woolf et al.| (1 200914 . Th e blue dashed line represents the fitting 
from IW»»iTetliTT(i2009l f . The blac k solid line represents our fit ¬ 
ting of G 22 4-58 B and data from IWoolf fc Wallersteinl ll200fil) : 
IWoolf et all (120091 1. The fit can be described with [Fe/H] = 
2.004 x CTiO/CaH — 1.894. Red dotted li nes shows t he bound aries 
of metallicity classes of M dwarfs from lLepine et all (120071 1. 

spectrum of A component we require the spectral data of a 
quality comparable with solar atlas. The under abundance 
found in oxygen may be explained as a real small number of 
atoms in the primary atmosphere due to its association in 
CO when the carbon abundance/presence is comparatively 
high in the atmosphere. Our picture could be better tested 
if were able to do more precise analysis of the CO bands at 
2.3 micron. 

The abundance inhomogeneities found can be most eas¬ 
ily explained in terms of the nucleosynthesis of halo stars at 
their formation epoch. The low iron abundan ce, i.e [Fe/H] 
= -1.9 2, together with large PM determined bv IZhang et al.l 
(120131) allows assignment of G 224-58 AB to the halo popu¬ 
lation. We detect in the atmospheres of G 224-58 AB, over¬ 
abundances with respect to the iron, not only of Ti, but 
also of other alpha-element elements. The study of similar 
systems, their physical characteristics as well as their mul¬ 
tiplicity, is of crucial importance for the investigation of 
binary and multiple system formation processes and their 
evolution through composition and time, including the for¬ 
mation of still hypothetical exoplanetary systems in the 
early epochs of Galaxy evolution. 
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